Understanding RHIC Collisions: Modified QCD 
fragmentation vs quark coalescence from a ther- 
malized fiowing medium 

Tom Trainor 

CENPA 354290, University of Washington, Seattle, Washington 98195, USA 
DOI: will be assigned 

The hydrodynamic (hydro) model applied to data from the relativistic heavy ion collider 
(RHIC) suggests that a dense QCD medium opaque to partons is formed in more-central 
Au-Au collisions. However, two-component spectrum analysis reveals a hard component, 
consistent with parton fragmentation described by pQCD which can masquerade as "radial 
flow." Minimum-bias angular correlations reveal that most scattered partons survive as 
"minijets" even in central Au-Au collisions. Such alternative methods quantitatively de- 
scribe spectrum and correlation structure via pQCD calculations. RHIC collisions appear 
to be dominated by parton scattering and fragmentation even in central Au-Au collisions. 

1 Introduction 

We wish to test the extent to which perturbative QCD (pQCD) can describe more-central A-A 
collisions at RHIC. Is a hydrodynamic (hydro) description necessary, or even allowed by data? 
Detailed arguments are presented in Refs. [TJ[HI31I1] with related material on hydro interpreta- 
tions in Refs. [51 IHl [Zl IB] ■ We describe pQCD fragment distributions (FDs) obtained by folding 
measured fragmentation functions (FFs) with a pQCD dihadron spectrum. We adopt a par- 
ton "energy-loss" model [9] to provide FD calculations which describe fragmentation evolution 
with A-A collision centrality. We also introduce a method to convert jet angular correlations 
to fragment yields and spectra. It is then possible to calculate the minijet contribution to the 
A-A final state and provide a comprehensive pQCD description of RHIC collisions. 

Measurement of azimuth quadrupole V2 via nongraphical numerical methods has led to 
strong conclusions about hydrodynamic evolution of a thermalized bulk medium with small 
viscosity. However, trends obtained from an alternative 2D V2 method suggest that the azimuth 
quadrupole is a novel QCD phenomenon carried by a small fraction of all final-state particles. 

2 Spectra and parton fragmentation 

The two-component model for pcr-participant-pair A-A spectra [ll llOj is 

_J_]_dnch ^ SNN{yt) + V HAA{yt,v) (1) 
npart yt ayt 

= SNN{yt) + vrAA{yt,v) HNN{yt), 
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Figure 1: First: Dijct spectrum for p-p collisions (solid curve) compared to an ab-initio pQCD 
theory result (bold dotted curve) and reconstructed jets (points) [5] , Second: Eq. ^ integrand 
for e+-e~ FFs, Third: Fragment distribution (solid curve) compared to p-p hard-component 
data (points) [10]. Fourth: Medium- modified FDs compared to Au-Au hard components 

where Snn is the soft component and Haa is the A- A hard component (with reference Hnn ^ 
Hpp). Ratio rAA = Haa/ Hjssm is an alternative to nuclear modification factor Raa- Centrality 
measure v = ^ntmary / npartidpant estimates the mean projectile-nucleon path length in A-A 
collisions. Soft component Snn by hypothesis remains unchanged with A-A centrality. For 
Glauber linear superposition of p-p (N-N) collisions (the GLS reference) spectrum hard com- 
ponent Haa HNNiut) also remains unchanged. In more-central A-A collisions HAA{yt,b) 
deviates from reference HNNivt), reflecting "medium modification" of parton fragmentation. 

We adopt the hypothesis that spectrum hard component Haa represents minimum-bias 
parton scattering and fragmentation in the form of a fragment distribution (FD). The pQCD 
convolution integral used to calculate fragment distributions is 

(P^h „ Tj f \ e(A7?)/2 [°" dadtjet 

— ^ 2TTytHAA[yt) ~ T / dymaxDrx{yt,yrnax)-i , (2) 

dyt drj cr„s„ A774^ Jo dymax 

where Dxx{yt, ymax) is the measured FF ensemble for collision system xx and dadijet/ dymax is 
the pQCD dijct (parton) spectrum [2]. The folding integral then predicts hadron spectrum hard 
component d^nt/dyt drj from parton pairs scattered into angular acceptance Ary. Comparisons 
with measured spectrum hard components in p-p and Au-Au collisions are shown in Fig. [T] The 
two-component model appears to provide a full description of A-A spectra in terms of pQCD. 

3 Correlations and parton fragmentation 

We want to determine the relation between minimum-bias jets (minijets) and the spectrum 
hard component (fragmentation). We integrate Eq. ([2]) over yt to obtain the hard-component 
yield 

2.HAA{b) « r ^...-e..(Ar;)/2 | r dymaxn.,,{ymaxM'^-p^] 0) 

I ^NSD ^m-K J L CTdijet Jo dymax J 

or 2-KHAA[b) = f{b)nchj{b), where the first factor is the pQCD jet frequency per NSD p-p 
collision [10] and the second factor is the minijet mean fragment multiplicity [llj . Jet angular 
correlations (same-side jet peak integrals in Fig. [5] - first two panels) can be represented by 
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Figure 2: Left panels: 2D angular correlations from 200 GeV Au-Au collisions for (resp.) more- 
peripheral and more-central collisions |12j . Third: Hard-component yield inferred from jet 
correlations (solid curve) compared to spectrum integrals (points). Fourth: Two-component 
total yield inferred from jet correlations (solid curve) compared to spectrum integrals (points). 

mean pair ratio j^(&) averaged over angular acceptance (A77, 27r) [TTJ[T2]. In terms of pQCD 
mean jet number nj{b) = niiin{b)Ar]f{b) and total multiplicity Uchib) per A-A collision in A?7 
we can write nch j{b) ~ nch{b) \/ p{b) /nj{b), which reveals the centrality dependence of the jet 
fragment multiplicity in A-A collisions. 

Figure [2] (third panel) shows spectrum hard component HAA{b) (solid curve) inferred from 
jet angular correlations via Eq. ([3]) . The open point is a spectrum estimate from Ref. [TU] (with 
Ar/ = 1). The solid points are derived from the "total hadrons" spectrum data in Fig. 15 (left 
panel) of Ref. [T]. Multiplying Eq. ([3]) through by v/2tt gives 



vHAA{b) is the hard component in the two-component spectrum model of Eq. IT}. Figure [2] 
(fourth panel) shows the two-component particle yield S^j^ + vHAAib) predicted by measured 
jet angular correlations (bold solid curve). Soft component s^vAf is by hypothesis fixed at 
~ 0.4 [2D density on (77, (j>)] for all A-A centralities. The solid points are the "total hadrons" 
data in Fig. 15 (left panel) of Ref. [T] divided by 27r. Minimum-bias angular correlations thus 
demonstrate that 1/3 of the final state in central Au-Au collision is contained in resolved jets. 

4 Spectrum structure and paradigm tests 

Jet quenching and Raa Spectrum ratio Raa presented as in Fig. |3] (left panels) suggests 
strong jet suppression and possible complete absorption of the majority of large-angle scattered 
partons in more-central Au-Au collisions at RHIC [13]. However, Raa presents a misleading 
picture because the ratio includes spectrum soft component S^n. In Fig. [3] (right panels) 
alternative hard- component ratio rAA reveals the true evolution of parton fragmentation with 
A-A centrality, demonstrating that suppression at larger pt is compensated by much larger 
enhancement at smaller pt. The number of jei-corre/aiec/ hadrons (hard-component multiplicity) 
in more-central A-A collisions increases dramatically compared to p-p collisions. 

Radial flow Radial flow is inferred by so-called blast- wave fits to pt spectra [M]. Such fits 
generally span a limited pt interval below 2 GeV/c attributed to "soft physics" described by 
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Figure 3: Left: Conventional spectrum ratios Raa for (rcsp.) pions and protons from 200 GeV 
Au-Au collisions. Right: Equivalent hard-component ratios rAA for (rcsp.) pions and protons. 



hydro models. But most of the parton fragments from jets fall below 2 GeV/c. It can be 
demonstrated that the spectrum evolution attributed to radial flow is quantitatively predicted 
by pQCD Inferred blast- wave parameters actually follow minijet systematica [3]. 



The baryon/meson ratio The baryon/meson (B/M) ratio in more-central A-A collisions 
is said to be anomalous in a pQCD context, suggesting that constituent-quark coalescence 
from a thermalized partonic medium may provide the dominant particle-production mechanism 
at intermediate pt [15| . However, examination of corresponding spectrum hard components 
suggests that parton fragmentation remains the fundamental hadronization process. 
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Figure 4: Left: Spectrum hard components for five centralities of 200 GeV Au-Au collisions 
for (resp.) pions and protons. Third: Soft and hard components for pions and protons from 
peripheral (i/ = 1) and central (v = 6) Au-Au collisions. Fourth: Proton/pion spectrum ratios. 



Figure |4] (left panels) shows spectrum hard components for identified pions and protons. 
The modes for both species in peripheral collisions are located at pt = 1 GeV/c reflecting a 
common underlying parton spectrum. For more-central collisions the pion mode moves down 
to ~ 0.5 GeV/c but the proton mode remains fixed on pt [1]. Both changes can be seen as 
FF modification as proposed in Ref. [9], but differently controlled in each case by the hadron 
fragment mass. Figure |4] (third panel) shows two-component parametrizations of pion and 
proton spectra accurate to about 10%. The parameterizations lead to proton/pion spectrum 
ratios (solid and dashed curves, fourth panel) which accurately describe the B/M "anomaly" 
data, details of which then correspond exactly to parton fragment distributions in the left 
panels. The structure in the fourth panel correspond exactly to that in Fig. [3] (right panels). 
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5 Correlation structure and paradigm tests 

Dihadron azimuth correlations Dihadroii azimuth correlations are intended to probe in- 
medium jet modifications in A-A collisions. A critical issue is subtraction of the combinatoric 
background modulated by an azimuth quadrupole measured by V2 and conventionally inter- 
preted as "elliptic flow." The absolute background offset is estimated by the ZYAM method, 
and V2 is obtained from published data derived from nongraphical numerical methods. ZYAM 
subtraction leads to two principal conclusions: (i) most partons are thcrnializcd in an opaque 
medium and (ii) the development of Mach cones by parton passage through the medium leads 
to distortion of the away-side azimuth peak (double humps) |16j . Close examination of the 
ZYAM subtraction method reveals that the "zero-yield-at-minimum" background estimate is 
not valid for the overlapping same-side and away-side jet peaks encountered in more-central 
A-A collisions. And published V2 data are typically overestimated due to strong jet (nonflow) 
contributions [Hl^. In effect, jet components are subtracted from other jet structure, and jet 
peak amplitudes are underestimated by an artifically high background offset estimate. 

Azimuth quadrupole The azimuth quadrupole conventionally interpreted as elliptic flow is 
measured by quantity V2 which can be defined in terms of nongraphical numerical methods |17| . 
Alternatively, V2{2D} can be obtained by fitting 2D angular correlations, and novel systematic 
trends then emerge [6]. Substantial differences between the two methods can be attributed 
to strong jet contributions to the former [7l|8j. The differences have implications for ZYAM 
subtraction and Mach-cone inferences [4]. The resulting energy and centrality trends of V2{2D} 
suggest that interpretation of the azimuth quadrupole as a hydro phenomenon is questionable. 




p,(lab) (GeV/c) y,(it,K,p) y,(it,K,p) y,(it,K,p) 



Figure 5: First: Published V2{pt) data for three hadron species. Second: The same data in a dif- 
ferent plotting format. Third: The same for As only. Fourth: Reconstructed quadrupole spectra 
for the three hadron species, illustrating a common source boost and emission parameters. 

In Fig.[S]the same published minimum-bias V2{pt) data for three hadron species are plotted in 
several formats. Plotted as V2{pt)/pt on transverse rapidity yt = ln[(pt -|-mt)/m/i] with hadron 
mass TO/i (middle panels) the V2 data reveal an apparent common source boost A?/to ~ 0-6 
and large deviations from a viscous- hydro theory based on Hubble expansion (R) [5]. When 
further multiplied by the single-particle spectrum poiut) the V2{pt) data (fourth panel) reveal a 
common quadrupole spectrum in the form of a Levy distribution on mt, much "colder" than the 
single- particle spectrum for most hadrons [71 [5] . Those results strongly suggest that "elliptic 
flow" in the form of the azimuth quadrupole is associated with a small fraction of the hadronic 
final state and results from QCD field-field interactions of low- a; gluon condensates [18]. 
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6 Summary 



Hydro analysis of RHIC data interprets partiele production below 2 GcV/c in terms of flow phe- 
nomena. The role of parton scattering and fragmentation is minimized. But model-independent 
analysis of spectrum and correlation structure reveals new fragmentation features quantitatively 
described by pQCD. pt spectrum hard components are manifestations of minimum-bias par- 
ton fragmentation quantitatively matched to minimum-bias jet angular correlations (minijets). 
Calculated pQCD fragment distributions accurately describe measured hard components. 

The reference for all fragmentation in nuclear collisions is the FD derived from measured in- 
vacuum e"'"-e~ FFs and the parton (dijet) spectrum for p-p collisions. Relative to the reference 
the spectrum hard component for p-p and peripheral Au-Au collisions is found to be strongly 
suppressed for smaller momenta. At a specific point on Au-Au centrality the hard component 
transitions to strong enhancement at smaller fragment momentum and suppression at larger 
momentum. 

Minimum-bias jet (minijet) correlations have been converted to absolute fragment yields 
which are found to comprise approximately one third of the final state in central 200 GeV 
Au-Au collisions, implying that almost all large-angle scattered partons down to 3 GeV parton 
energy survive as jet manifestations in the final state, albeit with some modification. pQCD 
calculations should be applied to all aspects of spectrum and correlation data to discover what 
is truly novel in RHIC collisions. pQCD describes almost all RHIC collision evolution — hydro 
interpretations are questionable. 

I greatly appreciate the hospitality of the ISMD organizing committee. This work was 
supported in part by the Office of Science of the U.S. DOE under grant DE-FG03-97ER41020. 
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